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We consider nuclear spin dynamics in a two-electron double dot system near the intersection of 
the electron spin singlet S and the lower energy component T+ of the spin triplet. The electron 
spin interacts with nuclear spins and is influenced by the spin-orbit coupling. Our approach is 
based on a quantum description of the electron spin in combination with the coherent semiclassical 
dynamics of nuclear spins. We consider single and double Landau-Zener passages across the S-T+ 
anticrossings. For linear sweeps, the electron dynamics is expressed in terms of parabolic cylinder 
functions. The dynamical nuclear polarization is described by two complex conjugate functions 
A* related to the integrals of the products of the singlet and triplet amplitudes CgCT + along the 
sweep. The real part P of A^ 1 is related to the S-T+ spin-transition probability, accumulates in the 
vicinity of the anticrossing, and for long linear passages coincides with the Landau-Zener probability 
Plz = 1 — e~ 27r "' ', where 7 is the Landau-Zener parameter. The imaginary part Q of A + is specific 
for the nuclear spin dynamics, accumulates during the whole sweep, and for 7 > 1 is typically 
an order of magnitude larger than P. P and Q also show critically different dependences on the 
shape and the duration of the sweep. Q has a profound effect on the nuclear spin dynamics, by 
(i) causing intensive shake-up processes among the nuclear spins and (ii) producing a high nuclear 
spin generation rate when the hyperfine and spin-orbit interactions are comparable in magnitude. 
Even in the absence of spin-orbit coupling, when the change in the the total angular momentum 
of nuclear spins is less than ft per single Landau-Zener passage, the change in the global nuclear 
configuration might be considerably larger due to the nuclear spin shake-ups. We find analytical 
expressions for the back-action of the nuclear reservoir represented via the change in the Overhauser 
fields the electron subsystem experiences. 

PACS numbers: 73.63.Kv,72.25.Pn,76.70.Fz 



I. INTRODUCTION 

Electron spin states in semiconductor quantum dots 
are investigated for their potential use as quantum bits 
in quantum computing architectures To this end, con- 
trol of the spin states and their couplings to the environ- 
ment is essential. In GaAs and InAs semiconductors, a 
major source of electron spin deco here nce is the coupling 
to the surrounding nuclear spins JaEHU Since the quan- 
tum dots are large compared to the interatomic spacing, 
each electron interacts with typically one million nuclei. 
Achieving control over this many-body interaction is a 
key for manipulating semiconductor quantum bits. 

In two electron double quantum dots, the singlet S and 
triplet To states define the elementary qubit. The cou- 
pling between these states is governed by the gradient 
in the longitudinal magnetic Zeeman splitting between 
the two dots. Controlling this coupling enables singlet- 
triplet qubit manipulations. Beyond the two-state S-Tq 
qubit operation, the gradient in the transverse magnetic 
Zeeman splitting between the two dots defines the cou- 
pling of the singlet S to the triplet T + and T_ states. 
Finally, the longitudinal magnetic Zeeman splitting de- 
termines the relative energies of the triplet states. This 
Zeeman splitting arises from the external field B and the 
nuclear spin background via the Overhauser field, and by 
changing the nuclear spin polarization the basic electron 
parameters can be tuned. 



Polarization of nuclear spins can be created and de- 
stroyed by flip-flop processes by pumping the elec- 
tronic states via time-dependent gate voltages. This 
has recently been investigated in many interesting 
experimental^^] anc j theoretical papers in double quan- 
tum dots in the regime of Pauli blockadeP^K' Experi- 
mentally, it has been demonstrated that an Overhauser 
field gradient of several hundred milli Tesla can be gener- 
ated and sustained.^ The dephasing time of the electron- 
spin qubits has been extended to more than 200 /xsP^ 
Because the dynamical interaction of an electron spin 
with a nuclear spin reservoir is enormously complicated, 
different theoretical efforts were focused on the various 
aspects of it. The two aspects most closely related to 
our paper are the theoretical modeling of the connec- 
tion between the generation of dynamic al nucl ear spin 
polarization at short and long time scales^^^ and the 
influence of the spin-or bit int eraction on the build-up of 
the nuclear polarization! 16 ! 17 ! 

The aim of this paper is to study in detail the elec- 
tron and nuclear spin dynamics as the system passes 
across a S-T + anticrossing. In GaAs and InAs quan- 
tum dots in an external magnetic field, T + is the lowest 
energy component of the electron triplet state because 
of the negative electron g-factor, g < 0. During a S- 
T + (or a T + -S) passage, electrons trade their spin with 
the nuclear reservoir, and multiple passages are used in 
creating a difference (" gradient" ) of the effective nuclear 
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(Overhauser) fields between two parts of the double dot 
that are used for qubit rotations. The study of a sin- 
gle passage (or two passages during a single cycle) pro- 
vides a firm basis for investigating events on longer time 
scales. Also, the progress in experimental techniques cur- 
rently allows, instead of averaging data over thousands 
of sweeps, to perform single-shot measurements and 
most recently such measurements have been achieved for 
double quantum dots.^ Also, the double dot dynamics 
during a single sweep manifests itself explicitly in beam 
splitter experiments!^ We expect the approach devel- 
oped in our paper to become a useful tool in discussing 
such types of experiments and, more widely, to facilitate 
better understanding and utilization of the nuclear spin 
environment in solid state based quantum computing. 

Specifically, we take into account the spatial distribu- 
tion of the hyperfine coupling between the electron and 
nuclear spins and compute the change in the topography 
of the nuclear spin polarization and the related changes 
in the gradient and average Overhauser fields governing 
the dynamics of the electron spin. These fields, that the 
electrons experience in the singlet and triplet states, de- 
pend on the spatial variation of the electron-nuclear cou- 
pling and we take this dependence into account. We em- 
ploy the Zener approach 21 and find analytically explicit 
expressions for the electron and nuclear spin dynamics 
during a single linear sweep and during cycles consisting 
of two linear sweeps. 

Let us give an overview of the main results. We express 
the whole electron and nuclear spin dynamics in terms of 
two complex conjugate functions A ± (Ti,Tf) depending 
on the initial and finite times (T,,T/) and the shape of 
the path between them. These functions are integrals 
of the products of the singlet and triplet amplitudes dur- 
ing the S-T+ passage. The real part P = RejA*} is 
the transition probability between the singlet S and the 
triplet T + states. The imaginary part Q = Im{A + } in- 
cludes basic information about the nuclear spin dynam- 
ics including the nuclear shake-ups. The Landau-Zener 
probability, Plz = 1 — e~ 2?r7 , where 7 is the Landau- 
Zener parameter, is the asymptotic value of P(T iy Tf) for 
a single sweep when Tj — > — 00 and Tf — > 00. Usually, 
all results are expressed in terms of Plz- Our approach 
provides a more detailed information about the nuclear 
spin dynamics away from the S — T + anticrossing. 

Oscillations of the transition probability P(Ti,Tf) as 
a function of its arguments reveal typical interference 
patterns. These oscillations are highly anharmonic for 
small Landau-Zener transition probabilities Plz ^ 1 
and might persist for a long time with a large ampli- 
tude for intermediate Landau-Zener transition probabili- 
ties Plz ~ 0.5. However, it is not typically the transition 
probability P that determines the nuclear spin dynam- 
ics. Instead, the other S-T+ quantity, Q is non less im- 
portant. While P is constrained to be in the interval 
< P < 1, there are no such constraints on Q and it 
is typically larger than P. We find that Q controls the 
shake-up processes among the nuclear spins. In the ab- 



sence of spin-orbit coupling, at most h of the angular 
momentum can be transferred to the nuclear spin bath. 
Given that there are around a million nuclear spins in the 
quantum dots, of which around a thousand are aligned 
initially, a change in one out of a thousand nuclear spins 
would have only a minor effect. However, the nuclear 
spins are allowed to interchange their spins during the 
S-T + passage without violating the conservation of the 
angular momentum. Although the interchange does not 
change the total nuclear spin angular momentum, the re- 
distribution of the nuclear spins can lead to considerable 
changes in the various gradient and average Overhauser 
fields that the electrons experience. This is because the 
Overhauser fields depend on weighted average values of 
the nuclear spin distribution with respect to the electron- 
nuclear couplings and not just the total nuclear spin. We 
find that such shake-ups are very sensitive to the initial 
nuclear spin distribution and that they are often much 
larger than the average nuclear spin production because 
Q is typically ten times larger than P. 

Furthermore, when the spin-orbit coupling competes 
with the hyperfine interaction and Q is considerably 
larger than P, then the Q-enhanced spin generation dom- 
inates for a generic direction of the nuclear spin polariza- 
tion and can become considerably larger than P. How- 
ever, after averaging over the direction of the transverse 
nuclear spin polarization, Q cancels and the results of 
Refs. |16I17| are recovered. 

Another finding is that even geometrically symmetric 
double quantum dots acquire asymmetric behavior be- 
cause of the spatial inhomogeneity of the hyperfine cou- 
pling. The sign of the asymmetry depends on £?, and its 
magnitude is largest close to the (0,2) or (2,0) configura- 
tion. The consequences of this B-controlled asymmetry 
for building nuclear field gradients are similar to that 
envisioned in Ref. 1151 for geometrically asymmetric dots. 

This paper is organized in the following way. In Sec. |Hj 
we describe the model of a double quantum dot that fol- 
lows the lines of Refs. [15122123] . We introduce the basic 
notations related to the electron-nuclear hyperfine inter- 
action and the nuclear dynamics induced by it in Sees. |III| 
and |IV[ respectively. In Sec. [Vj a linear Landau-Zener 
sweep is treated analytically and the time-dependence of 
the effective magnetic fields acting on the nuclei is dis- 
cussed in detail. Because Sec. [V] is rather technical, a 
reader interested in experimental applications can skip 
to Sec. VI where numerical data for the linear in time 
Landau-Zener sweeps and cycles are discussed. In Sec. 



|VII[ the back action of the nuclear spin dynamics on the 
Overhauser fields in the electron spin Hamiltonian is es- 
timated. Appendix [A| outlines the notations for electron 
spin operators. Appendix |B| discusses the spatial depen- 
dence of the hyperfine interaction. We demonstrate that 
even for two symmetric quantum dots, the hyperfine cou- 
pling acquires asymmetries controlled by the overlap in- 
tegral and the external magnetic field. Appendix [C] in- 
cludes two new identities for parabolic cylinder functions. 
We conclude and summarize our results in Sec. IVIIII 
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II. MODEL 

We consider two electrons in a double quantum dot. 
When the electron spin is conserved, the classification of 
the electron states as a singlet state S and three triplet 
(T„, v = 0,±1) states is exact. Spin-orbit interaction 
and the interaction with the nuclear spins mixes these 
states. We use the singlet and triplet stationary states 
as our basis. They are 



* s (l,2) =^(1,2)^(1,2), 
*t„(1,2) = *At(1,2)xt„(1,2). 



(la) 
(lb) 



where 1 and 2 denote the 1st and 2nd electron. The 
spin wave functions obey the symmetries xs(l>2) = 
— Xs(2, 1) as well as \t u (1, 2) = xt„ (2, 1) and are speci- 
fied in Appendix [B] The orbital wave functions -0s(lj2) 
and ^x(l>2) obey the symmetries 0s(l,2) = -05(2,1) 
and ^t(I) 2) = — 0t(2, 1), and we consider only the low- 
est energy orbital states so there are no additional quan- 
tum numbers labeling the orbital wave functions. 

The electrons interact with each other, external gate 
potentials, an external magnetic field, and with the nu- 
clear spins predominantly via the hyperfine interaction. 
The latter interaction, as well as spin-orbit coupling, in- 
duce transitions between the singlet and triplet states 
that we compute. The nuclei interact with the exter- 
nal magnetic field, the electrons through the hyperfine 
interaction, and with each other via the magnetic dipole- 
dipole interaction. The latter interaction affects the nu- 
clear spin dynamics on long time scales of around milli 
seconds, and we disregard it in what follows. However, we 
take into account (in a semiclassical Born-Oppenheimer 
approach and in the leading order in the large elec- 
tron Zeeman splitting) an indirect RKKY-like interac- 
tion between nuclear spins originating from the hyperfine 
electron-nuclear coupling (see Sec. VD| ). Near the ST + 
anticrossing it manifests itself at the scale of about 10 fia. 

Of central importance is the hyperfine electron-nuclear 
interaction 



4 /= A^^i(R r r ( )(i r sW), 



(2) 



where A is the electron-nuclear interaction strength, I 
numerates electrons and j nuclei, s(£) = \cr{t) are the 
electron spin operators in terms of the vector of Pauli ma- 
trices <t(£) for each electron £, and Ij are the nuclear spin 
operators. The electron and nuclear spin operators are 
dimensionless in our notations. Carets denote quantum 
mechanical operators and bold variables are vectors. 

In the 4x4 singlet and triplet space (S, T + ,Tq, and 
T_), the Hamiltonian that describes the electrons and 
their interaction with the nuclear spins can be written as 



H = 



v* e T ~fj-S 



(3) 



where the total electron spin S = s(l) + s(2). Addition- 
ally, the spin-orbit interaction induces terms in Eq. ^ 



that we discuss below. The nuclear spins are also af- 
fected by the external magnetic field through the nuclear 
Zeeman effect that we take into account below in the de- 
scription of their dynamics. However, we disregard the 
effect of the nuclear Zeeman energy on the equilibrium 
spin populations because of the high temperature of the 
nuclear spin bath. The eg and ct terms in the diagonal 
matrix elements of Eq. ^ describe the singlet and triplet 
energies in the absence of the nuclear and external mag- 
netic fields. They depend on the electrostatic gate po- 
tentials and the interactions between the electrons. The 
off-diagonal operator components = («+, — v^, — ?)~), 
are nuclear spin dependent (a superscript T denotes the 
transpose of a vector and the subscript n denotes that 
this coupling is due to the nuclear spins) 



(4) 



with a = (+, -, z), if = (i? ± iff) fy/2 are the trans- 
verse nuclear spin components, and the singlet-triplet 
electron-nuclear coupling coefficients 



(5) 



dependent on the positions R 3 of nuclei j. Roughly, pj 
varies from positive in one quantum dot to negative in 
the other. Therefore, and v z n represent differences 
in the effective nuclear magnetic fields in the two dots 
in the directions transverse and parallel to the external 
magnetic field, respectively. The effective splitting of the 
triplet states due to the external magnetic field B and 
the nuclei is —fj- S, where 



Vn 



(6) 



rjz is the electron Zeeman splitting in the field B || z, 
S is the spin-1 operator for the electrons (as defined in 
Appendix |A|) , and the position dependent coupling con- 
stants of the triplet states to the nuclei are 



dr#(r,R,)0T(r,R,). 



(7) 



This completes the description of the Hamiltonian that 
governs the coupling between the electron and nuclear 
spin dynamics. 

The ST + anticrossings arising due to v" and also 
the ST level splittings were investigated by th e beam- 
splitting technique^ and Rabi-oscillationj- 9 * 10 * 24 *, respec- 
tively. 



III. ELECTRON AND NUCLEAR SPIN 
DYNAMICS 

The Hamiltonian of Eq. ^ defines a many-body prob- 
lem of the coupled electron-nuclear dynamics. Our in- 
terest is in the dynamical nuclear polarization that is 
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achieved by changing the gate voltages in such a way 
that the electronic subsystem makes a transition from 
the singlet S to the lowest energy triplet T + state or vice 
versa. The many-body interaction can be simplified by 
employing the Born-Oppenheimer approach.^ The elec- 
trons are fast as compared to the nuclei. The electrons 
also interact with a large number of nuclei, around one 
million. These two features imply that the electron dy- 
namics is unaffected by the dynamics of a single nucleus 
and electrons see only a quasi-static configuration of all 
nuclei during a single ST + crossing. This motivates an 
ansatz where the wave function is separable into elec- 
tronic and nuclei partsP^ 

The electron dynamics can be solved from the Hamil- 
tonian of Eq. ([3| with the assumption that the nuclear 
spin operators can be approximated by their expectation 
values before the transition, v„ —¥ v„. The detuning 
energy e is defined as the difference between the triplet 
energy €t and the singlet energy eg, e = €t„ — eg, and 
is controlled by the variations in the gate voltages. We 
restrict ourselves to the limit of a rather large external 
magnetic field so that the splitting between the triplet 
states is larger than the magnitude of the off-diagonal 
matrix elements that mix the singlet and triplet states. 
When the separation between the energy levels is much 
larger than the matrix elements that mix the singlet and 
triplet states, the singlet and triplet states are well sepa- 
rated. The singlet-triplet matrix elements produce anti- 
crossings between the singlet and triplet levels when their 
energies are tuned to be close to resonance. Our focus is 
on situations where the system is tuned close to the S-T + 
transition as shown in Fig. [I] There, the energies of the 
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FIG. 1: Schematics of the singlet and triplet energy levels as 
a function of the detuning energy e = er — es close to the 
S-T+ anticrossing. The Zeeman splitting r\z — 1 is chosen as 
the energy unit, off-diagonal matrix elements are v ± = \v | = 
0.07. 

triplet states To and T_ are of the order the electron Zee- 
man splitting rf away from the energies of singlet S and 
triplet T + states, which is a large energy as compared to 
the S-T + anticrossing width. In this case, the electron 
dynamics can be approximated by the 2x2 dynamics for 
the singlet S and triplet T + amplitudes of the electron 
wave function. The reduction of the original 4x4 elec- 
tron dynamics problem to a 2 x 2 problem also facilitates 
finding an exact solution for the electron dynamics for 



linear sweeps and allows to reveal the role of the long 
time "tails" of the singlet and triplet amplitudes crucial 
for the nuclear spin dynamics. In the 2x2 basis, the elec- 
tron dynamics is described by the singlet cs and triplet 
ct + amplitudes that obey a Schrodinger equation 



H (ST + 



cs 

CT+ 



= ihdt 



cs 

CT+ 



with the Hamiltonian 



v~ 



(8) 



(9) 



where er + = £t — if and following Refs. |16|17j we have 
included the spin-orbit matrix elements v^ that couple S 
and T_|_ states into the total off-diagonal matrix elements 

v ± =vt + (10) 

While the coupling between S and T levels in GaAs 
double quantum dots is usually attributed to the hy- 
pcrfine interaction, spin-orbit coupling is inevitably 
present while difficult to evaluate quantitively for spe- 
cific devices.^ It manifests itself in spin relaxation ] 26 ! 27 ! 
level anticrossin gs in InAs single and d ouble dotspsHH 
and in the EDSfPH] both in GaAfpMland InAsP dou- 
ble dots. It is important to emphasize the existence of 
different mechanisms that couple the electron spin to the 
orbital degrees of freedom. They include the traditional 
(Thomas) spin-orbit interaction that couples the electron 
spin to the electron momentum and the Zeeman interac- 
tion in a inhomogeneous magnetic field B(r) that couples 
the electron spin to the elecron coordinate.^ In Ref. [32j 
the first mechanism dominated while in Refs. 36 and 1331 
different versions of the second one were important. We 
show in what follows that spin-orbit coupling also has a 
profound effect on the nuclear spin polarization produc- 
tion rate. 

By carrying out a unitary transformation of the orig- 
inal 4x4 Hamiltonian, it can be shown that the cor- 
rections to the reduced 2x2 Hamiltonian of Eq. Q are 
quadratic in the small ratio between v and the Zeeman 
splitting rjz provided the gate- voltage induced S-T + tran- 
sition is slow so that h(es — ^T)/Vz ^ 1- We assume that 
this criterion is satisfied. 

In turn, the dynamics of nuclear spins is driven by 
the effective magnetic fields Aj arising from the electron 
dynamics 



= Aj x !„ 



(11) 



where the components of the fields Aj acting on the nu- 
clei are the transverse = (A| ± «A^) /y/2 and longi- 
tudinal A| fields: 



A+ = Apjc s c^ + , 
At = A Pj c* s c T+ , 

&j = A (j\c T+ \ 2 ~ Vj(nZ), 



(12a) 
(12b) 
(12c) 
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and rjjinz) is the nuclear Zeeman splitting for the nu- 
cleus j. Because the dynamics of electron amplitudes 
(cs(t),CT + (t)) depends not only on the potentials on the 
gates but also on the nuclear spins through the matrix 
elements v , fields A.,- can be considered as dynamical 
RKKY fields. 

In the next section we show how the changes in the 
electronic states as they pass across the S-T+ anticrossing 
change the spatially dependent nuclear polarization. 



IV. DYNAMICAL NUCLEAR POLARIZATION 

We consider a situation where the changes in the gate 
voltages can induce a singlet S to triplet T + transition or 
vice versa, so that the total electron angular momentum 
may be increased or reduced by 1. In the absence of 
spin-orbit coupling, this implies that the change in the 
z-projection of the total nuclear spin equals the change 
in the elecron spin (but with the opposite sign). There 
is no conservation law for the spatial distribution of the 
nuclear spin. We are interested in how this change of 
angular momentum is distributed among the nuclei. As 
already mentioned above, the typical time scale for nuclei 
dynamics is long as compared to the time scale for the 
electron dynamics, in particular, with the singlet-triplet 
transition time. Let us denote the initial time of the 
sweep as T and the final time as Tf . We assume that the 
duration of the Landau-Zener sweep, Tf — T,, is short as 
compared to the typical nuclear spin precession time and 
take the nuclear dynamics into account as a perturbation. 
Also, since the total change of the angular momentum is 
of the order 1 , the typical change in the individual nuclear 
spins is much less than 1. With these assumptions, the 
change of a nuclear spin ALj = Ij(Tf) — Ij(Tj) during a 
Landau-Zener transition is 



(13) 



where the total effect of the electrons on the nuclei is the 



integrated effect of the magnetic splitting in Eqs. ( 12a 



(12b), and (12c 



r Tf fit 

T j (T f ,T i ) = J |A,(t) 



(14) 



In order to find explicit expressions for the dependence 
of the electron states on the effective field induced by the 
transverse nuclear spin polarization v^, it is convenient 
to make a transformation of the singlet and triplet am- 
plitudes 



c s = c s , c T+ = c T+ v /v ± . 
Then the Hamiltonian becomes real, and 



(15) 



where v± — \v |. Eq. (16) depends, in addition to the 
external magnetic field, on the absolute value of the com- 
bined effect of the nuclear spin induced transverse effec- 
tive field and spin-orbit interaction, but does not depend 
on its direction. 

In this basis, we can express the total effect of the (x, y) 
components of the effective field of Eq. ( 14 ) in terms of 



Tf =±iAp 3 A ± v ± /(2vl), 



(17) 



where the dimensionless functions A (Tj, Tf) are defined 

as 



A- = i2v ± j£* s {t)d T+ (t), 



(18) 



and A+ = (A-)*. This expression can be transformed 
by using the equation cV + = vj_ (ihdt — e s (t))c s following 
from Eq. (TBI), 



A- = -2 J ' dtc* s {t) d -^- -ilj' ^e s (t) \c s (t)\ 2 , 

(19) 

so that 



Re{A ± }=P = |c s (T; 



(20) 



is the transition probability P(T,Ty) from the singlet S 
to the triplet T + state. There is no such a simple relation 
between the imaginary parts of A ± and the transition 
probability, and this fact is important for the following 
discussion of the effect of the Landau-Zener sweeps on 
nuclei. However, we observe that when the Hamiltonian 



in the Schrodinger equation (16) is stationary, i.e., when 



the gate voltages are fixed and es and er + are constant 
in time, and the system is in an cigenstate of the Hamil- 
tonian of Eq. ( 16 ), the field c* s ct + is real implying a non- 
vanishing imaginary contribution to A*. The imaginary 
part of A* thus includes contributions that can be un- 
derstood in terms of RKKY-like static nuclear spin-spin 
interaction mediated by the electronic state, but this in- 
teraction also depends on the spin-orbit coupling. We 



will dia gona lize the stationary Hamiltonian of Eq. ( 16 ) 
in Sec. VD and relate the imaginary part of A* to the 
static electronic properties and show how this influences 
the dynamical nuclear dynamical properties. The imag- 
inary part of A 1 * 1 is central for the understanding of the 
dynamical nuclear polarization and we define 



Q = Im{A+} = -Im{A-}. 



(21) 



We also express the total effect of the field along z as 
rf = A( J A z /(2v ± ) - VKnZ)(T f - Ti)/h, (22) 

where 



cs 



= ihdt 



cs 

CT+ 



(16) 



A z = 2v ± 



Ti 



\ctM ■ 



(23) 
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Using Eqs. (12a I, ( 12b[ ) , and ( | 12c > , a s well as express- 



Eqs. (flTb and (23) in terms of A* 



mg J. J- and T z of 

and A z , we arrive at the spin production during a single 
S — > T+ transition both in the transverse 



An 



and the longitudinal components 



AI? 



A w + A+v + Pj P 



(24) 



(25) 



Next, substituting operators i n E g. Q by their semi- 
classical values v„ and using Eq. ( 10 ), we find the change 



in the z-component of the total nuclear spin, AI Z 
£AIJ, 



AP = -P+^ T [A-v-v+ + A+v+v-] , (26) 



AF = -J^( v -v++v+v-)-i^(v-v+-v+v-). (27) 



Note that the change in the z-component of the total 
nuclear spin is computed under the constraint that the 
transverse nuclear fields are before the sweep. 



Remarkably, AI Z of Eq. ( 26 1 only depends on the basic 
parameters of the Hamiltonian H ( ST +) of Eq. ^ and the 
shape of the sweep and does not depend on the detailed 
topography of nuclear spins. Therefore, the result is very 
general and convenient to use. In this respect, transfer 
of the longitudinal component of the angular momentum 
differs from the transfer of its transverse component that, 
according to Eq. (24 1, depends on the specific spin con- 
figuration. 



In the absence of spin-orbit interaction, 



f). the 



total change in the electron spin equals the transition 
probability P, as expected for a (partial) transition be- 
tween the singlet S and triplet T + states. Conservation of 
the z component of the angular momentum then dictates 
that the change in the z component of the total nuclear 
spin equals —P. Spin-orbit coupling breaks the conser- 
vation law for the angular momentum transfer from the 
electronic to the nuclear spin system since angular mo- 
mentum can be transferred to or from the lattice as well. 
Such processes manifest themselves in the second term 
in Eq. ( 26 ) . It depends on the relative phase between 



the spin-orbit and hyperfine interaction matrix elements. 
Furthermore, this term depends not only on the transi- 
tion probability P, but also on the imaginary parts of 
A*. Q acquires contributions not only from the part of 
the sweep near the anticrossing point but also from its 
long tails. As a result, the magnitude of Q can be much 
larger than P for certain classes of sweeps. This generic 



feature suggests that Q can be made large, and the spin- 
orbit coupling can strongly influence nuclear dynamics 
even when it is weaker than the hyperfine coupling. 

Our results confirm the prediction of Ref. [T7] that the 
spin-orbit coupling influences the nuclear spin generation 
rate profoundly. The quantity computed in Ref. [17] is the 
total change of the nuclear spin AI Z averaged over the 
phase of the transverse nuclear field v„ ■ This averaging 
annihilates the second term of Eq. (p7l) while the first 
term coincides with Eq. (9) in Ref. I17H-^ Since Q can be 
considerably larger than P, we expect enhancement of 
spin production rate in experiments performed at a fixed 
(while generic) values of vi~ . 



V. LINEAR SWEEPS - LANDAU-ZENER 
ELECTRON TRANSITIONS 

When the changes in the gate voltages are such that 
the difference in the energy betwen the singlet S and the 
triplet T + varies linearly in time, Eq. ( 16 ) reduces to the 



standard Landau-Zener problem. Because the Landau 
approach based on analytical continuation allows finding 
only the transition probabilities,^ we employ in the fol- 
lowing the Zener approach^ allowing finding explicit ex- 
pressions for the time dependence of electron wave func- 
tions that drives the coherent nuclear spin dynamics. We 
consider a transition from the singlet S state to the triplet 
T + state, but because of the symmetries of the Hamilto- 
nian the solution can also be used to find the wave func- 
tions that describe the transition from the triplet T + to 
the singlet S state. We derive this relation in Sec. |V C| 
Defining t = as the time when the energies e s and ct + 
of the singlet S and triplet T + are equal, we introduce 



/3 2 t/2h, e T = -/3 2 t/2h, 



(28) 



where /3 is a positive number with dimension of energy. 
This representation implies that the singlet state has the 
lowest energy at early (negative) times and the triplet 
state has the lowest energy for large final (positive) times. 
A natural time-scale is H/(3 so that Eq. (16 1 with t = 
tj3/h reads 



r/2 v/7 



cs 

C7V 



= id T 



cs 



where 



7 = {vx/py 



(29) 



(30) 



is the Landau-Zener parameter. When 7 is small, the 
transition probability from the singlet S to the triplet 
state T + is small. In the opposite limit, when 7 is large, 
the transition probability is close to 1. As above, we 
denote the initial time from where the sweep starts as Tj 
and the final time where it ends as T/. In dimcnsionless 
units, we have r,; = Ti/3/H and 17 = Tf/3/h. 

In order to determine the change in the nuclear spin 
polarization, we need to compute not only the transi- 
tion probability P, but also the singlet S and triplet T + 



7 



amplitudes, cs and ct + ■ Because the nuclear dynamics is 
controlled by the electron dynamics via the effective fields 



of Eqs. (12a), (|12b|), and (12c), explicit expressions for 



the amplitudes (cs(t),ct + (t)) should be found not only 
near the anticrossing point r = 0, but along the whole 
sweep, Tj < t < Tf. Therefore, it is necessary to em- 
ploy Zener's derivation of the Landau-Zener transition 
probability^ and complement it with a detailed infor- 
mation about the asymptotic behavior of the amplitudes 
and effective magnetic fields. 



Eliminating c s from Eq. ( 29 ) by substituting 



77(2+^ ' 



into its first row, we find 



d 2 T c T+ 



+ 7 



Then, by changing the variable r to 



(31) 



(32) 



(33) 



Eq. ( 32 ) transforms to 



d 2 z c T+ (z) + (n + i - \z^j ~c T+ (z) = 0, (34) 

where n = vy. This is the Weber equatio n 1 39 * 40 ! 
whose solutions are the parabolic cylinder (Weber) func- 
tions D n (z), D n (—z), D-i- n {— iz) and D-i- n (iz) of 
which only two are linearly independent. When ex- 
pressed as functions of the real argument r, they cor- 
respond to A 7 (e l37r/4 r), A 7 (-e l37r/4 r), £>_ 1 _ 47 (e l7r / 4 r) 



and D_i_i 7 (— e J7r / 4 r), respectively. In a similar way, we 
find the differential equation that the singlet amplitude 
obeys. Eliminating by substituting 



5T + = ^(-l+ W r) ZS, 



(35) 



into the second row of Eq. ( 29 ) and taking its complex 
conjugate, we find 



1 



(36) 



Hence c* s satisfies the same differential equation ( 32 ) as 



Op. ; its solutions are the Weber functions listed above. 



In 



Sec. |V A| we discuss the asymptotic behavior of the singlet 
S and triplet T + amplitudes that is critical for imposing 
the initial conditions and finding long time scale nuclear 
spin dynamics. 



A. Asymptotic Expansions 

For the following, the asymptotic behavior of the so- 
lutions in both limits, r — > ±00, is required. However, 
because the solutions appear in pairs, with opposite signs 
of r, it is sufficient to find their r > asymptotics. We 
note that the indeces of all above D-functions are imagi- 
nary or complex [ij or (— 1 — ry)] while the asymptotics of 
Refs. 1391401 are valid only for D n (z) functions with inte- 
ger indeces In what follows, we employ the asymptotic 
expressions from Mathematica 8 which are valid for arbi- 
trary complex indices. For large positive times r — > 00, 
they are 



'2it 



r(-*7) 

A 7 (-e i37r/4 r) « e *-7/Vr 2 /V 7 + 0(t~ 2 ), 

£>-l- i7 (e i7r/4 T) « e -W4 e7 r7/4 e -ir 2 /4 T -l- l7 + 0{t~ 3 ), 



-^/4 e -ir 2 /4 T -l-^ + 0( T - 2 ), 



V2tt 

r(i + i 7 ) ( 



-7r 7 /4 e ir 2 /4 r i 7 _|_ e j3ir/4 e -3ir 7 /4 e -ir 2 /4 r -l-i 7 _|_ Q^j--^ 



(37a) 

(37b) 
(37c) 

(37d) 



[asymptl-4] One can see that as r — > 00 the function 
_D_i_i 7 (e I7r7 / 4 r) vanishes as r _1 while the absolute val- 
ues of the three other D-functions saturate. We note 
that all asymptotic expressions for the D-functions in- 
clude two oscillatory factors. The Fresnel-type factors 
exp(±«T 2 /4) originate from the accumulation of the adi- 
abatic Schrodinger phases during a linear sweep, and the 
factors r ±n depending on 7 reflect the non-adiabaticity. 



It follows from Eq. ( |37c[ ) that for a sweep starting from 
the singlet S state at large negative initial time Tj, the 
function D_i_i 7 (— e" r ' r ' 4 T) should be chosen as one of 
the basis functions for the triplet T + state because it 
vanishes when r — > —00. We choose A 7 (e l37r7 ^ 4 r) as the 



second basis function. Then 



B. Infinite Limits and Asymptotics 



c T+ (r) = aVre-^C-i-vrC-e^/V) 

-— e-^DiJe^r), (38) 

where a and b are coefficients that depend on the initial 
time Ti. The overall phase factor as well as the factors yfiy 
and —1/^/7 have been chosen as a matter of convenience 
in the following transformation. One can check that b oc 
tT 2 for I tJ > 1. 



Eq. (361 implies that c* s , the complex conjugate of the 
singlet S amplitude, can be expressed in terms of the 
same Weber functions as the triplet amplitude ct + ■ An 
explicit connection between them can be found by em- 



ploying Eq. (31 1, and the expression for the singlet com- 



ponent c s can be further simplified b y usin g the standard 
recurrence relations for -D-functions! 39 * 40 ^ As applied to 
the D-functions of Eq. ( 38 ) , they read 



(I +id r ) A 7 (e i37r/4 r) = - 7 e^/ 4 D_ 1+J7 (e i3 ^ 4 r) 



(39) 



and 



(40) 

The D-functions of the right hand side of Eqs. (39) and 



( 40 1 differ from the D-functions of Eq. ( 38 1 , but are 



related to them by complex conjugation 

D_ 1+l7 (e^/ 4 r) = fz»_ 1 _ i7 (- e * r /V)l*, (41) 



D_ 47 (-e 



iw/4 



A 7 (e l37r/4 r) 



(42) 



Therefore, the general solution for the singlet ampli- 
tudes is 



cs (r) 



-(3 



+b 



^/ 8 A 7 (e l3 " /4 r) 



(43) 



As a consequence, the function A~ of Eq. (18) depend- 
ing on the product c* s (t)c.T + {t) and describing the re- 
sponse of nuclear spins to a Landau-Zener pulse can be 
expressed in terms of tw o fun ctions D_i_i 7 (— e l7r7 / 4 r) 

we consider the Landau- 



VB 



Zener scenario when the initial 
at ta - 



and D i7 (e l3 ^ 4 -r). In Sec._ 

electron state is prepared 
oo and the sweep runs to t/ — > oo, as well as 
the asymptotic behavior of effective fields c* s ct + at large 
but finite times |t| ^> 1. 



When the system is in the singlet state at early times, 
|cs(r — > — oo)| = 1 and ct + {t — > — oo ) = . then 6 = 
and \a\ 2 e*~i/ 2 = 1, as follow from Eq. (37b], and 



Ss(t) 



e *Vg-' r 7/4 



e- t3 ^ 8 A 7 (e j3 " /4 r) 



(44a) 



c T+ ( T ) = e ^g-7/4^ 



e- i37r/8 £>_ x _i 7 (-e i7r/4 r) 



(44b) 



where ^ is an arbitrary phase. For a finite but large initial 
time —Ti (ji > 0), this description remains satisfactory 
with the accuracy to the term s of the order t~ 2 in the 
singlet amplitude of Eq. (44a) and of the order t^ 1 in 
the triplet amplitude of Eq. ( 44b I . 



For completeness, let us also consider the situation 
when the system is in the triplet state T + at early times 
t — > — oo. Then it follows from Eqs. (37b) and (37c) that 
a = and e^/ 2 \b\ 2 /j = 1, so that 



cs(r) 

CT+ (t) 



^ [e-^/8 r;_ 1 _. 7 (_ e W4 T) ' 



_ _ e iip' e -7r 7 /4 



(45a) 

e- l3 " /8 A 7 (e j3 " /4 T)l , (45b) 



where <j)' is an arbitrary phase. 

We can now find the transition probability for the S 
T+ transition of Eq. Q. It is 



Plz = \c s (t -> -cx))| 2 - \c s (t -> oo)| 2 = 1 



-2-7T7 



(46) 

which is the celebrated Landau-Zener result. The trans- 
verse components of the effective field acting on the nu- 
clear spins are controlled by the product 



c s ct + = yffe-We-™*'* x (47) 



-7T7/2 -i3%/i 

A 7 (e i3 ' T/4 T)D_ 1 _ l7 (-e"/ 4 r) 



Its asymptotic behavoir following from Eqs. (37b) and 
d37dl is 



~c* s c t+ ^^+0{t~ 2 ) 



(48) 



for the early times r — > —00 and 



c s c T+ 



-i3ir/4_ 



sphx 
r(l + i 7 )' 



-3tt 7 /2 Jt 2 /2 2i~f 



T 2 ^ +0{T~ 2 ) 



(49) 
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for the late times r — > oo 
the second term of Eq. 



The absolute value of 
(49b is e^Vl - e- 2 ^ as 



can be checked by using the identity | T(l + ij) \ 2 — 
7Tj/ sinh(7T7). This result is easy to understand since it 
equals |cs||ct + | in the asymptotic regime r — > oo, where 



\or+\ 



1 — e 27r7 and \cs\ 2 — e 27r7 . The second term of 
Eq. ( 49 ) exhibits very fast Fresnel-like oscillations c %T / 2 
when t — > oo and does not contribute significantly to 
the integral A - of Eq. ( 18 ) describing the total effec- 



tive field applied to the nuclei as a result of the sweep. 
This factor originates from the accumulation of the phase 
exp {/ [es(t) — &r + (t)] dt/hj along the sweep. 

The origin of the coefficients in the 1/r terms in Eqs. 
(|48l and (49) can also be made quite transparent. By 



using the time-dependent Schrodinger equation (29), we 
find 



(id T + t) (c~s5 t+ ) = ^7 \c s (r)| 2 - \d T+ (r) 



(50) 



-oo, the amplit ude s 



Knowing that for early times, r — > 

1 and ct + = 0, we recover Eq. (48) 



approach |cg| 2 



For late times, |cg| — ct + - > — l + 2cxp(— 2ttj), which 
explains the 1/r term in Eq. (49). Furthermore, we note 
that in the leading order the operator (id T +r) annihilates 
the second term of Eq. (|49|). 
The integrals of Eqs 



(48) and (49) diverge logarith- 



mically when the integration limits approach ±oo. This 
means that while Plz of Eq. ( [46] ) and the total spin 
transfer AI Z of Eq. (26) (for v so ^ = 0) are controlled by 



the vicinity of the anticrossing point, the effective fields 
and shake up processes in the nuclear subsystem pro- 
duced by them are controlled by the global shape of the 
pulse. The same is true for AI Z when Vso^ 1 ^ 0. We note 
that while the presence of logarithmic terms is a general 
property of linear sweeps, they contribute to AI Z only in 
the presence of spin-orbit coupling. 



C. Reverse sweep from the triplet T+ to the singlet 

S. 

Let us relate the reverse sweep, starting in a triplet 
state T + and sweeping to a singlet state S, to the S —> T + 
sweep elaborated above. Since now the rates of the 
change of the singlet S and triplet T + energies have 
the signs opposite to the signs in Eq. (28), the dynam- 



ical equations for the amplitudes (cs,c~t + ) differ from 
Eq. (16) by the interchange c~s *H> 5x+. Furthermore, 
for a T + — > S transition, the system was initially in the 
triplet T + state, hence, the singlet S amplitude vanishes 
at the early time. Therefore, the initial conditions are 
also cs ■<-> ct+ interchanged as compared to the S — >• T + 
sweep. This implies that their product transforms as 
c* s ct + — > (c* s ct + ) , and A* — > — (A*) according to 
Eq. ( [i~8| }-. In other words the transition probability 
P = Re{A } changes sign, but the imaginary parts 
Q = Im{A + } remain unchanged. The change of the 



sign of Re{A ± } is obvious because of the S O T + in- 
terchange, so that the longitudinal component of the an- 
gular momentum transfer changes sign. However, the 
effective field ImjA 1 * 1 } does not change, and this indi- 
cates that the imaginary components of A* should add 
during a S — » T + — > S cycle. 

In conclusion of this section, for linear sweeps the di- 
mensionless function A~(Ti, Tf) that reflects the effect of 
a single Landau-Zener sweep on nuclei diverges logarith- 
mically when Tj — > — oo and Tf — >■ oo. In Sec. |VI[ we 
discuss in more detail the dependence of A^ 1 (Tj ,Tf) on 
the limits (Tj,Ty) and the Landau-Zener parameter 7. 



D. Adiabatic Regime 

Some more insight on the long-r tails of the pr odu cts 
c* s ct + comes from the stationary solution of Eq. (16)P^ 
For a large detuning 8 = £t + — es from the S — i+ an- 
ticrossing, when |r| 3> 1, the stationary solution of Eq. 
( |16| provides an adiabatic approximation to the singlet 
and triplet amplitudes. Note that we still assume the du- 
ration of the sweep is short as compared to the nuclear 
Larmor precession time. 

Then the eigenenergies of the electronic states of the 
Hamiltonian of Eq. ( 16 ) are 



1 



05/2)' 



(51) 



and at the lower branch of the energy spectrum the prod- 
uct of the amplitudes equals 



c s c T+ 



vj2 



(52) 



Here the oscillatory r-dependent phase factors cancel be- 
tause cs and ct + belong to the same eigenvalue. It im- 
mediately allows calculating the transverse components 
= AgjC~sCT + v± / v ± an d Aj = AgjC* s CT + v~ /v± of A 
and the effective fields from Eqs. (12a) and (12b) 



transverse components vanish as v±/5 w 



j 

The 



len j d j /v± 



00. Similarly, the longitudinal component found from 



Eq. ( 12c I equals 



AQ 
2 



5/2 



(S/2Y 



Vj(nZ)- 



(53) 



Far from the intersection, when d/vj_ — > —00 and the 
eigenstate is almost a pure triplet, A| — > —AQ — rjj^ nZ y 
In the opposite limit, when S/v± — > 00 and the eigenstate 
is almost a pure singlet, A| — > —r]j( n z)- The point 5 = 
has been identified as "spin funnel" in Ref. 1321 

In the adiabatic limit, the fields Aj acquire the usual 
meaning of RKKY fields with a nuclear dynamic time 
scale of t ~ h/Aj. Near the level anticrossing point 8 = 0, 
Aj ~ Ano /N where no is the concentration of nuclei and 



10 



N is the number of nuclei in the dot. With Ano « 10 
cV and N « 10 6 , t « 10/zs. 

For a slow linear sweep between = 
<5 — > /3r, one finds from Eqs 
which, according to Eq. 



Tf and t f, with 



18 I and ( 52 I the quantity 
21 1, result in 



<3(o) = 4 7 In 




(54) 



and from Eq. (20) we find Pr a \ = 0. The results for P( a ) 
and Q( a ) hold with logarithmic accuracy; the subscript 
(a) indicates that they were derived in the adiabatic ap- 
proximation. In the same way, one can check that c* s cq\ 
of Eq. ( 52 ) is in agreement with the 1/r terms of Eqs. 
and |49J). 

Applying Eq. (52) to a nonlinear dependence S = 6(t), 
one easily concludes that A ± converges if 5(t) is super- 
linear and diverges by some power law if it is sublinear. 

Equation ( 52 ) implies important consequences for the 
nuclear spin dynamics under the condition of time- 



independent detuning. Indeed, it follows from Eqs. (10), 
( 11 ) - ( 12b I, and ( 52 ) that the rate of change of the total 



nuclear spin is 



dl z 



Therefore, time-independent detuning results in produc- 
ing a magnetization I z that increases linearly in time as 
long as the parameters of the electronic Hamiltonian re- 
main unchanged. This generation of spin magnetization 
by time-independent electrical bias is possible because 
the time-inversion symmetry is violated by a strong ex- 
ternal field B producing Zeeman splitting of the electron 
triplet state, and the simultaneous presence of hypcrfine 
and spin-orbit interactions. The magnitude of the ef- 
fect reaches its maximum at 5 = 0, when the system is 
brought to the center of the ST+ anticrossing. The time 
scales of the parameter change can be estimated similarly 
to Sec. | VHP | Under the usual conditions, the shortest 
of them corresponds to the precession of in the ex- 
ternal field. These conclusions seem to agree with the 
observations of Rcf. 031 



VI. S -> T+ SWEEPS AND ROUND CYCLES 



Complex functions A ± (T i ,Tf) of Eq. (18) describe 
the effect of a sweep on the nuclear spins. As seen 
from Eqs. (20) and ( [26| ), the probability of the electron 
S —> T + transition P is completely controlled by the 
real part of A ± , P — Re{A ± }, while the angular mo- 
mentum transfered to the nuclear system AI Z depends 
both on the real and imaginary parts of A*. Imaginary 
parts of A^ are always present but manifest themselves 
in the nuclear spin accumulation only when there are two 
competing mechanisms of the electron spin transfer, hy- 
pcrfine and spin-orbit. 



In this section, we first present data on the dependence 
of A* on the integration limits and the Landau-Zener pa- 
rameter 7 obtained by numerical integration of Eq. ( 29 1 , 
and then develop an analytical approach for describing 
the oscillatory dependence of the transition probability 
P on the cycle length. 



A. Linear sweeps 

We begin with linear S — > T + sweeps of Sec. W\ For 
such sweeps, we denote the initial time — Tj (tj > 0) and 
the final time Tf (r/ > 0) so that the duration of the 
sweep is + Tf. To reduce the number of parameters, 
we assume Tj = Tf. Transition probabilities P(t/) are 
plotted in Fig. 2(a) as a function of the sweep half-time 
Tf for two values of 7. While for large 77 both curves sat- 
urate to the Landau-Zener probabilities Plz of Eq. (46 1, 
oscillations of P{jf) are very pronounced. 



They decay 

at a rather long time scale, and their shape cannot be 
described by a single characteristic time. We attribute 
the oscillations to the interference pattern between two 
spectrum branches and estimate their period t osc from 
the Schrodinger exponent exp(— iv±t/K) in the anticross- 
ing point, what results in r osc 7 _1 / 2 . The rate of their 

(55) 

decay is controlled by the passage time Hv± //3 2 across the 

1/2 



avoided crossing that results in a decay time Td ec ~ 7 
Finally, we arrive at a rough estimate of the transient 

}. Actually, this only is 



regime T tr ~ max{7 1 / 2 , 7~ 1 / 2 
a lower bound on r tr . The saturation takes a longer 
time and the difference in the shapes of the 7 = 1 and 
7 = 0.1 curves deserves more comments. The 7 = 0.1 
curve strongly resembles plots of Fresnel integrals, and 
we attribute the oscillatons to the e lT I 2 factors in the 
asymptotics of Eq. (49). With increasing 7, the pat- 



terns of oscillations are getting less regular due to the 
second oscillatory factor r 2 * 7 in the asymptotics of c*ct + • 
The switching of regimes happens at 27r7 w 1 as is seen 
from the expression e~ 27r7 for the Landau-Zener transi- 
tion probability. 

In agreement with the asymptotics found in Sec. |VB[ 
the imaginary parts of A displayed in Fig. 2(b,c) exhibit 
a behavior quite different from the behavior of their real 
parts P. They increase nearly logarithmically with 77, 
with weak oscillations superimposed on this monotonic 
growth. Their magnitudes increase with 7, and for 7 rs 1 
and Tf ~ 10 they are by one order of magnitude larger 
than P. Therefore, even with a moderate spin-orbit cou- 
pling, the imaginary parts of A ± are expected to con- 
tribute essentially to the spin transfer AI Z of Eq. ( p6"| . 
This contribution should not only change the magnitude 
of AI Z but also smoothen its 77-dependence. 

In Fig. |2(b)[ we also plot Q = Im{A+} for 7 = 1 by 



using the approximate adiabatic expression of Eq. (54) 



to compare it to the exact numerical result. Apart from 
some details of the behavior for early and late times, 
which are expected, we see that the dominant contribu- 
tion to Q can be explained in terms of the adiabatic field 
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P=Re[A ± ] 
1.0 h 




Tf(=Td 



t/(=t0 



Q=Im[A + ] 




T/(=Ti) 



(c) 



FIG. 2: (a) Transition probability P = Re{A ± (r / )} for a 
linear sweep starting in the S state at the initial time — n 
and ending at the final time tj = Ti plotted as a function 
of the half-sweep time r/ for two values of the Landau-Zener 
parameter 7. Full (blue) line 7 = 1, dashed (red) line 7 = 
0.1. The anticrossing point is passed in the middle of the 
sweep at time r = 0. The full (blue) lines in (b) and (c) are 
Q — Im{A + (r/)} for 7 = 1 and 7 = 0.1, respectively. In (b) 
and (c), the dashed (green) lines are the adiabatic solutions 
of Eq. JHil. 



of Eq. (54). Fig. |2(c) provides a similar comparison, 
but for a faster sweep with 7 = 0.1. Even in this situa- 
tion, the adiabatic approximation is a reasonable starting 
point for describing the basic shape of Q of Eq. [2l] 

The above analysis of linear sweeps, together with 
the arguments of Sec. |VD allow to make some conclu- 



sions about the generic (non-linear) S-T+ sweeps as well. 
Imagine the sweeps with the rate unchanged near the 
anticrossing but increasing away from it. As long as the 
speed-up happens at times r > r tr (this inequality should 
be fulfilled strong enough), the probability P — RejA*} 
changes only modestly, while the long time tails of the 
products c* s (t)ct + {t) contributing to Q — Im{A + } are 
cut-off. Thus, increasing the sweep rate away from the 
anticrossing reduces Q and might have a profound effect 
on AI Z . However, its specific magnitude depends on the 
values of a number of parameters such as , vf a , r tr , and 
the speed-up time. 



B. Cyclic linear sweeps 

Round sweeps are of the most practical interest for 
experiment, and their detailed shapes are nontrivial be- 
cause of the oscillating tails of RejA*} of Fig. 2(a). 
Therefore, we provide below the data on for two dif- 
ferent round sweeps starting in the singlet states S at 

n < 0. 

Fig. [3] presents data for a round sweep of the total 
duration of 4r/ that includes the sweep of Fig. 2 from 
Ti = —Tf to 77 and the backward sweep that begins im- 
mediately after the end of the forward sweep. According 
to Eq. ( |20| , P = Re{A ± } displays the probability of 
S->T + transition. Remarkably, Fig. 3(a) shows that for 
7 = 1 the decay of P is rather long and includes deep 
and irregular oscillations. For 7 = 0.1, P(t/) shows a 
wide maximum at 77 ~ 2, and the following oscillations 
without any visible decay up to Tf = 10. In this case, 
a double dot in the linear sweep regime resembles a res- 
onator of a length decreasing as tJ 1 . We expect that 
first peaks can be resolved experimentally, e.g., in beam 
splitter experiments^ while higher peaks should merge 
into a background with P w 0.5. Using first sharp peaks 
for ultrafast spin operation is highly tempting. 

As distinct from P = RelA*}, Q = Im{A+} of 
Fig. 3(b) is a nearly monotonic function of 77 for 7 = 1 
(with irregular oscillations superimposed), and is about 
10 for Tf — 10. Therefore, it can heavily contribute to 
AI Z . However, ImjA 111 } is small and strongly oscillates 
at 7 = 0.1. 

To demonstrate the effect of the tunneling process near 
the anticrossing point, in Fig. [4] are plotted the data for 
a cycle that begins in the S state at — 77, reaches the 
anticrossing at r = 0, and then runs immediately back 
with the same speed until 77 with r» = 77 . Compari- 
son of Figs. 3(a) and 4(a) for 7 = 1 shows quite similar 
patterns of the oscillations of P{Tf) that are more regu- 
lar in Fig. 4(a). However, the patterns for 7 = 0.1 are 
rather different demonstrating essential decrease in the 
spin transfer. The magnitudes of Q — Im{A + } are small 
in both cases, but their Tf dependences are rather differ- 
ent. 
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P=Re[A ± ] 
l.Oh 




JmMJ Tf (= T j) 



Q=Im[A + ] 




t/(=t ; ) 



RelA*] 

l.Or 




2 4 
-Im[A-](=Im[A + ]) 



Tf(=-Ti) 




Tf(=-Td 



FIG. 3: (a) Transition probability of a S — > T+ transition 
P — Re{A ± } and (b) the imaginary part Q = Im{A + (r/)} 
for a round sweep plotted versus 77 (one fourth of the sweep 
time). The first part of the sweep is the same as the sweep of 
Fig. [2j and the second part sweeps in the opposite direction 
with the same speed immediately after reaching the turning 
point. Full (blue) lines 7 = 1, dashed (red) lines 7 = 0.1. 



C. Analytical theory of the probability oscillations 

We can explain the oscillations in the transition prob- 
ability as a function of the total duration of the cycle 
employing the analytical results in Sec. |V| The forward 
sweep from — to the turning point r m gives rise to the 
singlet and triplet amplitudes of Eq. (44 1. Assuming 



T m ^ 1 and Tj 3> 1, and employing (37a I and (37d), the 
singlet and triplet amplitudes at the turning point T m are 



(56) 



FIG. 4: Real part P (a) and imaginary part Q (b) of the 
function A + (r/) for a round sweep plotted vs half-sweep time 
Tf. First part of the sweep, starting in the S state at n < 
r < 0, stops in the anticrossing point at r = and runs 
immediately in the opposite direction until tj — \n\. Full 
(blue) lines 7 = 1, dashed (red) lines 7 = 0.1. 



Sec. VC 



for the backward sweep (eg ^ , ) 



the dynamical equations for the amplit udes 
differ from Eq 
by the interchange cs < — > Ct_, ■ In order to ni 



(16) 
ate 



contact with our results in Sec. [Vj we change the 
time r — > r — 2r m for the backward sweep. Using the 
interchange eg < — > c^ + , it follows from Eq. (44b) that 
for the triplet T + channel the ratio of the final and initial 

amplitudes along the backward sweep is Cg + /a^f^ — 
^ye-i^/BD^ (_ e "/4 r/ ) / [ e -i3^/8 A7 (.gW*/^)]^ 



T + in the superscript of c^ T+ 



indicates the channel. In 



the limit r/ 3> 1, Eqs. (37d) and (37c) imply that this 
ratio equals 



£ (fS) 



e -i37r/4 e -7r 7 /2 



r(i + *7) 



(57) 



Here the superscripts indicate that we started in the sin- 
glet S state and carried out a forward linear sweep. The 
phase of the early time singlet state is arbitrary and is 
omitted because it only modifies the overall phase of the 
wave function and does not influence the final result for 



the probability. The amplitudes of Eqs. ( 56 1 and ( 57 1 are 



derived under the assumption that eg = 1 and cV + = 
at time r = — t<. 

Next, we consider the backward sweep and include 
the contributions from two channels passing through the 
T + and S states at the turning point. As discussed in 



~(bT + ) ,JfS) 



1%: 



e -i37r/4 e - 



■7T7/2 



J/4 il i-i 



r(i + »7) 

Similarly, by using the interchange cs < — > ct + , it follows 
from Eq. (45b) that for the singlet S channel the ratio 



of the final and initial amplitudes along the backward 



~(bS) ,~(fS) 

sweep c s /Cg 



D 



&7 



3itt/4 



Tf) /Dr 



a 347T/4 
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the limit Tf ^> 1, Eqs. (37d) and (I37ch imply that the 
ratio of the singlet amplitudes after the backward sweep 
equals 



-(bS) ,-tfS) 



>/2g > «e^^e l(T /" T "^ /4 T) 7 T- l T. (59) 
The singlet amplitude at the final time r/ after the cycle 
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of duration (r» + r m ) + (r m + r/) is a sum of the contribu- 
tions coming from both channels, Cg 0t ^ = Cg T+ ^ + Cg S \ 
Finally, 



C S 



(1 - Plz) + P L ze li>(Tm) 



(60) 



where the Landau-Zener transition probability Plz for a 
single-passage is defined by Eq. (46 ) and the phase ^(r m ) 
at the turning point r TO is defined as 



'T 4i7 e" 



M / 2 r(-i 7 )/r(z 7 ). (6i) 

The dependence of the transition probability P = 1 — 
|cs(t/)| on the position r m of the turning point is 

P(r m ) = 4F LZ (1 - P LZ ) sin 2 0(r m )/2, (62) 

where $(r m ) is the Stuckelberg phaseP^^It is acquired 
between the two passages and includes both the adiabatic 
and non-adiabatic (7-dependent) parts. From Eq. (62) 
we can make several observations that are consistent 
with the numerical data of Fig. [3j First, when r/ >• 1 
and t, ^> 1, P does not depend on the initial and final 
times. The transition probability only depends on the 



Landau-Zener probability Plz of Eq. (46 1 and the turn- 
ing point T m . This means that the oscillations of P(r m ) 
are a robust feature of a coherent double passage across 
a Landau-Zener anticrossing. The transition probability 
oscillates around the average value 



Pav = 2P L z(l-PLz). 



(63) 



For fast sweeps Plz < 1 so that P oscillates between 
and APlz- For slow sweeps Plz is close to 1 and the 
probability oscillates between and 4(1 — Plz)- The 
maximum in the oscillation amplitudes is achieved at 
Plz = 1/2. When P LZ = 1 - e~ 2 ^ = 1/2 (7 « 0.11), 
the transition probability P oscillates between and 1. 
The amplitudes of the oscillations are smaller for all other 
values of 7. This is exacly the behavoir we see in the nu- 
merical plots. One more remarkable feature of Fig |3][a) , 
that all oscillations pass through P — 0, is also reflected 



by Eq. (621 



Oscillatory patterns of 7 = 0.1 curves in Figs. 2(a) and 
3(a) show strikingly different behavior. In Fig. 2(a), the 
amplitude of oscillations decreases with t/ , and P gradu- 
ally approaches its Landau-Zener limit Plz- On the con- 
trary, in Fig. 3(a) the oscillations, after some transitional 
period, acquire a stationary amplitude. Eqs. (62 1 and 



(63) clarify the origin of this behavior typical of double 
pass ages across the anticrossing.SSfUl Indeed, P av of Eq. 



(2) 

(|63|) is a Landau-Zener probability P L % f° r a double pas- 



sage across the anticrossing that can be derived directly 
by the above two-channel procedure with quantum am- 
plitudes substituted by probabilities, see Ref. [35] There- 
fore, suppression of these long-time scale oscillations and 



(2) 

approaching the double-passage Landau-Zener limit P L % 
are only achieved when the decoherence is taken into ac- 
count, and can allow measuring decoherence times. 

In conclusion, prolonged oscillations of the electronic 
amplitudes (cs,ct + ) are a generic property of the co- 
herent electron dynamics during the single- and double- 
passages across the S-T + anticrossing. Their amplitudes 
and durations are controlled by the Landau-Zener param- 
eter 7 and by dephasing on longer time scales, and the 
patterns are rather different for the single- and double- 
passages. 



VII. BACK ACTION OF NUCLEAR SPIN 
DYNAMICS ON OVERHAUSER FIELDS 



The Hamiltonian H of Eq. (|3j) describing the electron 
states depends on the Overhauser fields created by the 
spatially dependent nuclear spin configuration. The elec- 
trons experience the nuclear fields v" and rj n of Eqs. Q 
and ([6]), where the first represents the components of the 
effective difference magnetic field in the dots, and the sec- 
ond represents the induced average magnetic field. When 
going through the S-T+ transition, the electrons will ex- 
perience a change of these nuclear Overhauser fields. It 
is a unique property of Eq. (26) for the change in the 
total longitudinal nuclear spin AI Z that it expresses a 
global property of a double dot in terms of the parame- 
ters of the electronic Hamiltonian and does not depend 
of a specific configuration of nuclear spins. For differ- 
ent elements in the Hamiltonian H, we calculate their 
mean-square values as well as their variances. 

The expression for the change of the total z compo- 
nent of the nuclear spin of Eq. (27) makes the role of 



Q explicit due to the mediation of spin-orbit coupling. 



With 



ii. the total spin transfer is protected by 



the momentum conservation law and Q manifests itself 
through shake-up processes in the nuclear spin reservoir 
respecting the conservation of the total angular momen- 
tum. The electron dynamics induces changes in the nu- 
clear spin configuration that in turn induce changes in 
the in the diagonal and off-diagonal elements of the elec- 
tron Hamiltonian (|3j). In what follows, we compute these 
changes. 



A. Changes in Overhauser fields 

Electrons experience an effective Zeeman splitting in 
the Overhauser field of r)j of Eq. (16]). The associated 
change in the z-component of fjj, Ai)^ = —AJ2j CjAiJ, 



Ar)* 



A*_ 
2v\ 



(64) 



In the multicycle regime^the field of Eq. (64) has been 
measured by Petta et al.™ and by Foletti et alP by the 
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shift in the position of the ST+ anticrossing. In contrast 
to AI Z , the change Afj^ in the longitudinal field depends 
on the detailed nuclear spin configuration and on the spa- 
tially dependent electron-nuclear couplings pj of Eq. ^ 
and Q of Eq. Q. 

The singlet-triplet terms and in the Hamiltonian 
H of Eq. |3| are sums over all nuclear spins. ST level 
splittings characterized by v„ were measured in Ref. |2H 
and a number of follow-up papers, and ST + splittings 
described by t)„ in Ref. [2Ql The changes in these terms 
during a cycle are Aw" = A Y,j Pjlf- By using Eq. (24 1, 
we find changes in the components a = ± that couple - ! S 
to T± 



A 2 
2v ± 



P.I- 



3 

Vj(nZ) 



(65) 



and, by using Eq. (25), in the component a = z coupling 
S to T 



Av z = - 



A 2 
2v]_ 



(66) 

We note that while only produces a longitudinal Over- 
hauser field mixing S and To, Av^ includes operators 
and therefore mixes S and T + belonging to our 2x2 
subspace. 

In the next sub sections, mean values and variances of 
these operatores are computed. 



B. Constraints and mean values 

While nuclear spins are distributed in the bath ran- 
domly, the magnetization fluctuations controlling 
electron dynamics during the cycle impose on their values 
the constraints 



(67) 



adding also a constraint related to w*. To simplify cal- 
culations, we consider below the nuclear spins Ij as ran- 
dom Gaussian variables that are normalized, in the ab- 
sent of constraints, as (Ijlji ) = + l)§jj'8\\>, with 
A = (x, y, z). Then the mean values of Ij are 

3 Ilffdl^HipWvk-A^pjd},) 

(68) 

where V(Ij) are Gaussian probabilities, (v^,v^) are de- 
fined as = (v„ ± u^)/v2, and the denominator se- 
cures the normalization of the probabilities under the 



Using the integral representation for <5-functions 
1 



6(x) 



2tt 



(69) 



(70) 



multiple Gaussian integrations of Eq. ( 68 ) result in 

{If) = Pjvt/{AR 2 ), (/?) = Pj v*J(AR 2 ), 

where R n = Y]j p™ are determined by the spatial depen- 
dence of the electron-nuclear coupling constants. Sub- 
stituting these expressions into Eqs. (64) and (65), we 



arrive at the corrections to the nuclear field experienced 
by the electron spin during the sweep. 



(A<) = -APAR'JR 2 , 

_ P 2 (j, ai 
S and T + components 



(71) 



where R' 3 = J2j PjCj, an d the Overhauser field mixing its 



(A<) = -AFAR 3 /R 2 , 



(72) 



with AI Z of Eq. (26) 



We see that both the changes in the longitudinal dif- 
ference field Avi*. and the longitudinal average field Ar/ Z 
are proportional to the change in the total nuclear spin 
AI Z . It follows from Eqs. (Ts]) and ^ that pj typically 
have opposite signs in both dots while Q > everywhere, 
hence, R' 3 > 0. Therefore, with A > 0, the sign of (Ary^) 
(the change in the mean Overhauser field building in the 
double dot) is opposite to the sign of AI Z , in agreement 
with Eq. Q. The sign of (Av^) is defined by the sign R 3 
that depends on the choice of electronic basis functions 
(see Appendix B), therefore, it is not uniquely defined 
with respect to AI Z . 

The magnitudes of Ai] z and Av^ are of the order of 
AI Z Ano/N per cycle, i.e., about AI Z /\^N of the mean 
values of ?] z and v z t . For vf a — 0, AI Z = —P, hence, 
I AI Z |< 1. However, it is seen from Figs. 2(b) and 
3(b) that Q is an order of magnitude larger than P when 
7 > 1. Therefore, when v so ^ 0, the conditional ex- 
pectation values (Ar^) and (Av^) should experience Q- 
enhancement through the Q-enhancement of AI Z , and 
i] z and can change by about 1% per cycle. 

The mean values of the transverse components of v„, 
calculated in a similar way from Eq. (65), are 



(Au?) = A 



± 2 



2v\ 



R ^A±± l A V ^ R *A z 
i? 2 2v ± R 2 



where 



V(nZ) 



is a mean value of 



Ti)/n 

Vj(nZ) 



(73) 

over all nuclear 



species. Because different species are distributed ran- 
domly at the scale of atomic spacings, they self-average 
in the linear approximation over Tf — Tj, and we accept 
that all of them have the same absolute values of the an- 
gular momenta, Ij = I. While the first term is compara- 
ble in the magnitude to Eq. (72 ), the two last term might 



constraints of Eq. (67) 



be much larger because they increase with the sweep du- 
ration. However, Eq. ( 73 1 includes changes both in the 
amplitude and the phase of Av^, and the latter might 
not be essential when solving Eq. ( 16 ) that only depends 
on v±. We come back to this term in Sec. |VII D| 
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C. ST+-pulses induced interdot shake-ups 

Let us explain the importance of the variance in the 
spin production by considering the total nuclear spins in 
the left and right dots. Average values of different oper- 
ators calculated in Sec. IVH Bl were based on the condi- 
tional mean values (I") of nuclear spins I" of the order 
of N' 1 ' 2 that are small compared with their root mean- 
square values. Therefore, calculating the mean-square 
values of all operators and their variances is important 
for estimating the widths of statistical distributions. 

We begin with the differences in the spin polarizations 
of the left and right dots, L and R, that are critical for 
spin manipulation. While division of a double dot into 
its left and right parts holds only when the overlap in- 
tegral is small enough, cf. Appendix B, the results are 
instructive. Splitting Eq. Q into sums over L and R, we 
define partial sums 



V nL(R) — A 



jeL(R) 



(74) 



Their sums are v% and are a subject to constrains of 
Eq. (67). However, their differences 



J nR 



(75) 



are free of any constraints. Using Eq. (25), the change in 
the left-right polarization difference is 



ai Ir = -^r( A ~ w ~< + A+v+ u-)- (76) 



When averaged over an unpolarized spin reservoir, its 
mean value vanishes, (A/£^) = 0, and the mean-square 
value equals 



((A/; 



LR) 



A 2 n Q 
6t>^ 



1(1+1) | A | 2 yV(R)cf 3 R, (77) 



with p(R) of Eq. ^ and 

| A \ 2 =P 2 + Q 2 . 



(78) 



A simple estimate of the right hand side of Eq. ( 77 ) re- 
sults in | A | 2 . Therefore, the asymmetry of spin pumping 
of the left and right dots is Q-enhanced whenever Q ^> P, 
in particular, when v so — and P < 1. We attribute this 
enhancement to shake-up processes resulting in multiple 
spin flips per each "pure" injected nuclear spin. These 
processes are random, and it is not clear for now h ow 
they influence inhomogeneous spin distributions 1121111 

The detailed spatial patterns of spin generation at long 
time scales are a subtle subject and are related to the 
spatial variation of the electron- nuclear couplings p(Rj) 
and C(R-j') calculated in Appendix |B| With mean values 



of Ij of Eq. ( 70 1, spatial distribution of AI Z is related to 
AI Z as Alf = {p 2 , 

„2 



(p 2 /R 2 )AI z . The left-right asymmetry 
m Pj originates either from the geometric asymmetry of 



the double dolfl^ or from the L-i?-overlap of the electron 
density, cf. Appendix B, and produces a regular differ- 
ence in the I z generation rate. While the results depend 
on the specific distribution of nuclear spins and the <S*- 
Tq mixing,^ the mechanism of Q-enhancement is quite 
general whenever 7 > 1. 



D. Mean-square values and variances 



Mean values of Sec. |VIIB| were evaluated over an unpo- 
larized nuclear spin bath and estimate the mean rates of 
the change of the different parameters. However, the esti- 
mate of the shake-up rate of Sec. | VII C] demonstrates that 
calculating variances of these random variables can pro- 
vide additional, and sometimes even more valuable, in- 
formation about the magnitudes of the expected changes 
during a cycle. The conditional probability distributions 
are so wide that the mean value is not very representa- 
tive. In this section, we evaluate variances of the basic 
nuclear fields. 

We begin with calculating the mean-square values. Be- 
cause all nuclear fields of Eqs. ( 64 ) - ( 66 1 are linear in the 



momenta I" , mean values of the quadratic forms in them 
include integrals that differ from Eq. (68) by substitut- 



ing Ij~ either by (I^) 2 or by Ijl^, with j ^ j'. While the 
latter terms are smaller in the parameter 1/N <C 1, they 
have a higher statistical weight. Summing all terms, one 
arrives at length expressions for ((A^) 2 ) and ((Au^) 2 ) 
that we do not present here. Instead, using the mean 
values of Eqs. (71) and (72), we present the variances 



defined as Var{|| = (<f) - (£) 



.4 



Var{A^} =| A | 2 ^/(I + l)[i^ - (R> 3 ) 2 /R 2 ], (79) 
where R' A = J2APjQ) 2 , an d 



Var{A<} =| A | 2 ^-I(I + l)[R 4 - (R 3 ) 2 /R 2 ]. (80) 



(71) and (72) shows Q 



Comparison with Eqs. 
enhancement even when v so = (hence, when AI = 
— P), the effect that manifested itself already in Eq. (77). 
This means that the nuclear spins with I z away from the 
mean conditional expectation values (Ip respond to the 
sweeps stronger than the spins with H — (Ij)- Also, 
this enhanced sensitivity is due to the spatial distribu- 
tion of pj and Cj because with p 7 = const and (j=const 
the brackets in Eqs. (79) and (80) vanish. By the or- 
der of magnitude, both quantities experience changes of 
about AAno/N per cycle; with A sa 10 and N ps 10 6 , this 
suggests changes about 1% per cycle. In other words, 
around 10 spins interchange their directions during one 
passage. 

Calculating (A(v^[v^)) results in a simple equation 



(81) 
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because the contributions of the two last terms of Eq. 



( 65 1 cancel. In absence of spin-orbit coupling, this im- 
mediately suggests (A(vj_)) — PAv*R 3 /R 2 . Under these 



conditions, large terms in Eq. ( 65 ) reflect only the change 
in the ph ase of ii„ that does not influence dynamical 



equations (16), and the relative change in (v^) 2 = v£v~ 
is only about N~ x / 2 . 

However, in presence of spin-orbit coupling the dy- 
namics of spin amplitudes (cs,ct) is controlled by 
rather then ■ Mean value of A(w^), calculated by us- 
ing Eqs. (f65l) and (Tfoj) , is 



(A(v + v-)) = (A(«+«-)) 
^^-v-v++A + v+v-] 



+ i{ v n v to- v t v so) 



A R' 



¥< 3 >- 3i >-si:is A, 



(82) 



where first term is defined by Eq. (81). Physically, sec- 
ond and third terms in Eq. ( 82 ) take into account the 



angle between i>+ and v~^ a in the complex plane, and are 
proportional to the product v±v^. With v± ~ v^, rela- 
tive corrections coming from the second term are of the 
order A/y/N per cycle. The third term is usually much 
larger because it increases linearly with the pulse dura- 
tion AT = Tf — Tj. It includes two contributions of which 
first is due to the Zeeman precession of nuclei and sec- 
ond due to the Knight field and is proportional to the 
integral of | ct+ | 2 - While the magnitude of the second 
contribution depends on the shape of the pulse, the ratio 
of these terms is roughly fj^ n gy/(Ano/N) and they be- 
come comparable at B ~ 1 mT. This indicates that first 
contribution to the third term usually dominates. With 
v± ~ vj~ and fj( n B) ~ 10 mT, the Zeeman term results in 
(A(v + v~)) ~ 0A(v + v~) for a 0.1/iS linear sweep. This 
is much larger th an t he correction to the same quantity 
estimated in Eq. (81) and to ((Ai^) 2 ) having the same 
scale. The effect in InAs should be much larger than in 
GaAs because of the stronger spin-orbit coupling. 

The above estimates indicate that, because of the 



terms in Eq. ( 65 ) linear in the pulse duration, spin-orbit 



corrections to transverse matrix elements are essentially 
larger than the corrections to the longitudinal ones. 

In Eq. ( 82 ) , Zeeman precession of nuclei manifests itself 
in ((v + v~)) only through spin-orbit coupling. The effect 
is much stronger when estimated through the variance of 



and we estimate it for 



= when 



,± - „± 



Disregarding two first terms in Eq. ( 65 ) , the calculations 



similar to those performed above when deriving Eqs. ( 79 ) 
and (80) result in 



Var{A(«>-)} 



6 



x {v+v-^MiTf-TA/h] 2 , (83) 



where rj 2 n ^ is the mean-square value of rjj( n B) ■ It follows 



is the nuclear spin precession with a characteristic time 
of about a microsecond at B ~ 10 mT. It is about two to 
three orders of magnitude shorter than the corresponding 
time for estimated above. 

It is also instructive to compare this estim ate with a 



much longer time for V- following from Eq. (81). The 



latter estimate was found with the nuclear configuration 
of Eq. ( 70 ) that reflects the mean- values of nuclear spins 
under the constraints of Eq. (67). In a narrow region 



of the phase space around these mean values dy nam ics 
of is strongly suppressed. The estimate of Eq. (83) is 



much more representative because it represents the entire 



phase space compatible with the constraints of Eq. ( 67 1 . 



A similar type of the behaviour of was discussed above 
as applied to Eq. ( 80 ) . 



from Eq. ( 83 1 , the dominating mechanism of changing 



VIII. CONCLUSIONS 

We have studied the dynamics of the electron and nu- 
clear spins near ST + avoided crossings in double quan- 
tum dots. While adopting the traditional approach based 
on the hierarchy of time scales, with a slow nuclear and 
fast electron dynamics, we employed a quantum descrip- 
tion of the electron spin and coherent dynamics of nu- 
clear spins, and investigated the time-resolved patterns 
of single and double Landau-Zener passages through the 
anticrossing point. They are described by two complex 
conjugate functions A ± depending on the initial and fi- 
nite times (Ti,Tf) and the trajectory of the sweep, with 
A~ proportional to the integral of the product c* s (t)cT(t) 
of the complex amplitudes of the S and T+ states. Their 
real parts P = R^A 1 * 1 } are proportional to the S-T+- 
transition probability and for one-side sweeps oscillate 
at small time scales when the system is close to the 
anticrossing and saturate at long time scales. For lin- 
ear sweeps, we find the singlet and triplet amplitudes 
in terms of Weber /^-functions (parabolic cylinder func- 
tions); the long-time asymptotic limit of P equals the 
Landau-Zener probability Plz = 1 — e~ 27r7 . For round 
trips, the system also experiences long-term Stiickelberg 
oscillations. The first sharp oscillations might be uti- 
lized for ultrafast electron spin operation, while the de- 
cay of the oscillations can provide information about de- 
phasing rates. It is important that the imaginary part 
Q = Im{A+} that acquires contributions from the elec- 
tronic states at a wide time scale and accumulates with 
time (it diverges logarithmically for linear sweeps) has 
a profound effect on the dynamics of the nuclear spins. 
When the Landau-Zener parameter 7 > 1, Q is typically 
one order of magnitude larger than P. Therefore, in pres- 
ence of the spin-orbit coupling violating the angular mo- 
mentum conservation, Q may become the major factor 
controlling the angular momentum transfer to nuclei. In 
particular, this mechanism is efficient for excursions in- 
cluding a stay near the anticrossing point. Generically, 
A = (P 2 + Q 2 ) 1 ! 2 controls the shake-up processes that 
exchange angular momentum between the left and right 
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dots. With Q ^> P, it is Q that plays a dominating role 
in these angular-momentum exchange processes. Because 
the mechanism that plagues many experimental efforts of 
building considerable polarization gradients remains un- 
known, it is a challenging question whether and how the 
shake-up processes contribute to it; unfortunately, only a 
theory including multiple passages can resolve it. We also 
estimated changes in the Overhauser fields during a sin- 
gle cycle and concluded that the transverse components 
are more volatile than the longitudinal ones. 

We are grateful to B. I. Halperin, C. M. Marcus, I. 
Neder, and M. Rudner for stimulating discussions and 
comments on the manuscript. E. I. R. was funded by 
IARPA through the Army Research Office, by NSF under 
Grant No. DMR-0908070, and in part by Rutherford 
Professorship (Loughborough, UK). 



Appendix A: Spin Operator 

We use the following convention for the spin-1 operator 



s z = 




These operators satisfy the 



(Al) 
(A2) 
(A3) 

commutation 



relations 



Si , Sj 



= ieijkSk, where e^k is the Levi- 



Civita tensor, as well as S% + Sy + <§| = 2. 



Appendix B: Simple Model 



The singlet part of the spin wave function is 

Xs(l,2) = ^(|ti)U 3 >-Ui)|t 3 » (Bl) 

and the three triplet components of the spin wave func- 
tion are 



XT + (l,2) = |tl)|t2>, 



(B2a) 



XTo (l,2) = -=(|ti)|| 2 ) + |li)|t2», (B2b) 



V2 

Xt_(1,2) = Ui)|! 2 ). 



(B2c) 



We will in this section discuss the spatial dependence 
of the hyperfine coupling constants pj of Eq. ^ and Q of 
Eq. ([7]). In a simple model, the electron wave functions 



near the S-T+ anticrossing are 

[V>L(l)Vii(2) + Vz,(2)Vfl(l)], (B3) 



^5(1,2) = 0031^(1)^(2) 

sin v , 



V2 



^r(l,2) = -^[^(l)V«(2) 



Vl(2)Va(1)], (B4) 



where L denotes the left and R the right dot, and the 
angle v depend on the Zccman energy rjz- The normal- 
ization coefficients in ( B3 1 and ( B4 1 are exact under the 



assumption that the functions ipx, and i/'fl are orthonor- 
malized. 

Let us illustrate the spatial dependence of the electron- 
nuclear coupling constants p of Eq. ^ and £ of Eq. 
([7]) for a simple model of a quantum dot. We assume 
the electrons are in the lowest orbital harmonic oscillator 
state. The Cartesian coordinates, the wave function is 
ip(x,y) — exp [— (x 2 + y 2 )/l 2 ] / {ly/2/ir), where I is the 
size of each quantum dot. We have two quantum dots 
that are separated at a distance d, one at x = —d/2 and 
y = and the other at x = d/2 and y = 0. We form an 
orthonormal basis set based on the functions tp(x—d/2, y) 
and ip(x + d/2, y). In this basis, we compute p(x, y) and 

C(x,y)- 

We plot in Fig.[5]the electron-nuclear couplings p(x, y) 
and C( x yy) f° r y — as a function of x when v = 0.1 
and v = 7r/2 — 0.1. The spatial distribution of the 
singlet-triplet coupling p(x,y) depends on the angle v. 
When v is close to n/2, there is a nearly equal prob- 
ability for electrons to be located in the left and right 
dot for both the singlet and triplet states. Then the 
singlet-triplet coupling p(x,y) is nearly antisymmetric 
around x = 0, p(x, y) ps — p{— x, y ) [th e sign of p(x, y) de- 
pends on the sign choice in Eq. (B4)]. When v is small, 
the electrons are in the singlet state (0, 2) in the right 
dot, so that p(x, y) passes through zero inside the right 
dot (for x > 0). Therefore, even for two symmetrically 
shaped dots, the S-T + electron-nuclear coupling can be- 
come asymmetric because of the overlap of the left and 
the right dot wave functions. The asymmetry depends 
on v controlled by the external magnetic field. 

The triplet-triplet electron-nuclear coupling C( x ,y) 
does not depend on v and is a symmetric function of 
x for the two symmetric quantum dots. 



Appendix C: Two identities for the parabolic 
cylinder D-functions 



Using the solution of Eq. ( 44 ) for ct + (t) and c$ (r) and 



the normalization conditionfes (t) | + |&r + (t)| = 1, we 
arrive at an identity 

7 I £U_ i7 (-e i7r / 4 r) | 2 + I A 7 (e i37r/3 r) | 2 = e^' 2 (CI) 

relating absolute values of two D-functions at arbitrary 
real values of r and 7. 
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FIG. 5: The spatial variation of the electron-nuclear couplings 
p{x,y = 0) (a) and ((x,y — 0) (b). In (a), the red (full) 
curve is for v = 7r/2 — 0.1 and the blue (dashed) curve is 
for v = 0.1. The size of the dots is I — 50 nm and the 
separation between the dots is d = 100 nm. The overlap 
integral between the left and the right oscillator wave function 
is 0.1. It is the most striking feature that the overlap between 
the wave functions induces asymmetry of the p(x, y) even in 
geometrically symmetric double dots. The asymmetry reaches 
its maximum when the system is close to the (0, 2) state. 



Next, it follows from Eq. (291 that 



d T (| c s | 2 - | ~c T+ | 2 ) = -2*^/7 (csCT + - c s c* T+ ) . (C2) 



Integrating it over r and using Eqs. (44) and (46), we 
find 



/oo 
dr Im{e- i3 "/ 4 J D^ 1 _ ! ; 7 (-e"/ 4 r)A 7 (e 43 ' r/4 T)} 
OO 

(C3) 



sinh7T7 e _ 7r7/2 
7 



The integral of the real part of the integrand diverges. 



While we could not find these identities for D n (z) func- 
tions with complex (imaginary) indeces n and the argu- 
ments directed along diagonals in the complex z planes 
in any of mathematical sources, we checked them numer- 
ically. 



1 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Rev. Mod. Phys. 70, 1217 
(2007). 

2 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

3 J. Levy, Phys. Rev. Lett. 89, 147902 (2002). 

4 I. A. Merkulov, Al. L. Efros, and M. Rosen, Phys. Rev. B 
65, 205309 (2002). 

5 A. Khaetskii, D. Loss, and L. Glazman, Phys. Rev. B 67, 
195329 (2003). 

6 S. I. Erlingsson and Yu. V. Nazarov, Phys. Rev. B 70, 
205327 (2004). 

7 W. Yao, R.-B. Liu, and L. J. Sham, Phys. Rev. B 74, 
195301 (2006). 

8 L. Cywinski, W. M. Witzel, and S. Das Sarma, Phys. Rev. 
Lett. 102, 057601 (2009). 

9 S. Foletti, H. Bluhm, D. Mahalu, V. Umansky, and A. 
Yacoby, Nat. Phys. 5, 903 (2009). 



10 H. Bluhm, S. Foletti, D. Mahalu, V. Umansky, and A. 
Yacoby, Phys. Rev. Lett. 105, 216803 (2010). 

11 H. Bluhm, S. Foletti, I. Neder, M. Rudner, D. Mahalu, V. 
Umansky, and A. Yacoby, Nat. Phys. 7, 109 (2011). 

12 C. Barthel, D. J. Reilly, C. M. Marcus, M. P. Hanson, and 
A. C. Gossard, Phys. Rev. Lett. 103, 160503 (2009). 

13 G. Ramon and X. Hu, Phys. Rev. B 75, 161301 (R) (2007). 

14 M. Stopa, J. J. Krich, and A. Yacoby, Phys. Rev. B 81, 
041304(R) (2010). 

15 M. Gullans, J. J. Krich, J. M. Taylor, H. Bluhm, B. I. 
Halperin, C. M. Marcus, M. Stopa, A. Yacoby, and M. D. 
Lukin, Phys. Rev. Lett. 104, 226807 (2010). 

16 M. S. Rudner and L. S. Levitov, Phys. Rev. B 82, 155418 
(2010). 

17 M. S. Rudner, I. Neder, L. S. Levitov, and B. I. Halperin, 
Phys. Rev. B 82, 041311 (2010). 

18 H. Ribeiro and G. Burkard, Phys. Rev. Lett. 102, 216802 
(2009). 



19 



19 W. Lu et al., Nature (London) 423, 422 (2003; J. M. Elzer- 
man et al, Nature (London) 430, 431 (2004); T. Meunier 
et al., Phys. Rev. B 74, 195303 (2006); S. Amasha et al., 
Phys. Rev. Lett. 100, 046803 (2008). 

20 J. R. Petta, H. Lu, and A. C. Gossard, Science 327, 669 
(2010). 

21 C. Zener, Proc. Roy. Soc. (London), Series A, 137, 696 
(1932). 

22 J. M. Taylor, H. A. Engel, W. Dur, A. Yacoby, C. M. 
Marcus, P. Zoller, and M. D. Lukin, Nat. Phys. 1, 177 
(2005). 

23 J. M. Taylor, J. R. Petta, A. C. Johnson, A. Yacoby, C. 
M. Marcus, and M. D. Lukin, Phys. Rev. B 76, 035315 
(2007). 

24 J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 
Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and 
A. C. Gossard, Science 309, 2180 (2005). 

25 M.P. Now ak, B. Szafran, F .M. Peeters, B. Partoens, and 
W. Pasek, |arXiv:1102. 10021 

26 D. V. Bulaev and D. Loss, Phys. Rev. B 71, 205324 (2005) 

27 A. V. Khaetskii and Y. V. Nazarov, Phys. Rev. B 61, 12639 
(2000). 

28 A. Pfund, I. Shorubalko, K. Ensslin, and R. Leturcq, Phys. 
Rev. B 76, 161308 (2007). 

29 S. Takahashi, R. S. Deacon, K. Yoshida, A. Oiwa, K. Shi- 
bata, K. Hirakawa, Y. Tokura, and S. Tarucha, Phys. Rev. 
Lett. 104, 246801 (2010). 

30 E. I. Rashba and V. I. Sheka, in Landau Level Spectroscopy 
1991 (North- Holland, Amsterdam), p. 131. 

31 V. N. Golovach, M. Borhani, and D. Loss, Phys. Rev. B 
74, 165319 (2006). 

32 K. C. Nowack, F. H. L. Koppens, Y. V. Nazarov, and L. 



M. K. Vandersypen, Science 318, 1430 (2007). 

33 M. Pioro-Ladriere, T. Obata, Y. Tokura, Y.-S. Shin, T. 
Kubo, K. Yoshida, T. Taniyama, and S. Tarucha, Nature 
Physics 4, 776 (2008). 

34 S. Nadj-Perge, S. M. Frolov, E. P. A. M. Bakkers, and L. 
P. Kouwenhoven, Nature, 468, 1084 (2010). 

35 E. I. Rashba, J. Supercond. Novel Magnetism, 18, 137 
(2005). 

36 E. A. Laird, C. Barthel, E. I. Rashba, C. M. Marcus, M. P. 
Hanson, and A. C. Gossard, Phys. Rev. Lett. 99, 246601 
(2007). 

37 We are grateful to M. Rudner for his help in establishing 
this connection. 

38 L. D. Landau and E. M. Lifshitz, Quantum Mechanics 
(Pergamom, London, 1965) §90. 

39 E. T. Whittaker and G. N. Watson, A Course of Modern 
Analysis (New York, MacMilland, 1945). 

40 Higher Transcendental Functions, ed. by Erdelyi (Mc. 
Graw-Hill, New York, 1953), Chapter 8. 

41 In Zener's paper^, the aymptotics of RefP^ have been 
used. Fortunately, this did not influence his final result 
for the Landau-Zener transition probability. 

42 J. R. Petta, J. M. Taylor, A. C. Johnson, A. Yacoby, M. D. 
Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Phys. Rev. Lett. 100, 067601 (2008). 

43 S. Foletti, J. M artin, M. Dolev, D. Mahalu, V. Umansky, 
and A. Yacoby, |arXiv:0801.36T3l 

44 E. C. G. Stiickelberg, Helv. Phys. Acta 5, 369 (1932). 

45 E. Shimshoni and Y. Gefen, Ann. Phys. (NY) 210, 16 
(1991). 

46 S. N. Shevchenko, S. Ashhab, F. Nori, Phys. Rep. 492, 1 
(2010). 



